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Abstract 

We study renormalization group flows in strongly interacting field theories with 
fermions that correspond to transitions between a theory without Lorentz invariance 
at high energies down to a theory with approximate Lorentz symmetry in the infrared. 
Holographic description of the strong coupling is used. The emphasis is made on 
emergence of chiral fermions in the low-energy theory. 


1 Introduction 

Recent years have witnessed a revival of the idea [1, 2] that Lorentz invariance (LI), instead 
of being a fundamental property of Nature, may emerge only as an approximate symme¬ 
try of the low-energy physics. An important motivation to explore this possibility comes 
from the proposal by P. Hofava [3] that abandoning LI improves the ultraviolet behavior of 
gravity. In this proposal, departure from LI at large energies allows to postulate anisotropic 
(Lifshitz) scaling [4] which renders the gravity theory power-counting renormalizable. There 
is a version of Hofava-Lifshitz gravity [5, 6] which is viable and passes the existing experi¬ 
mental constraints [7]. On the phenomenological side, breaking of LI has been invoked in the 
construction of interesting long-distance modihcations of gravity [8, 9, 10]. These provide a 
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playground for the analysis of possible deviations from general relativity and may be relevant 
for cosmology [11, 12]. 

Emergence of an effective Lorentz symmetry at low energies is a fairly common phe¬ 
nomenon in condensed matter physics, see e.g. [13, 14, 15, 16, 17, 18]. Also, it is not hard to 
construct particle physics models with this property. Indeed, the original works [1, 2] consid¬ 
ered Yang-Mills theory and quantum electrodynamics with an addition of Lorentz violating 
terms of dimension four^ in the Lagrangian and showed that the coefficients of those terms 
vanish along the renormalization group (RG) flow towards the infrared (IR). This result 
persists in the Lorentz violating extension of the full Standard Model [19]. Furthermore, a 
general argument [20, 21] shows that LI hxed points of RG flows are IR stable with respect 
to GPT-invariant deformations and thus have non-empty basin of attraction in the space 
of GPT-invariant non-relativistic theories. All theories belonging to this basin of attraction 
exhibit emergent LR. In two space-time dimensions it has been shown [22] that under broad 
assumptions interacting hxed points with isotropic scaling are necessarily LI. 

However, this scenario faces an important challenge when confronted with the precision 
tests of LI in particle physics experiments [23, 24, 25]. In a weakly coupled theory the 
running of the coefficients in front of marginal Lorentz violating operators is logarithmic. 
If they are of order one^ at a certain ultraviolet (UV) scale A*, say, of order Planck mass, 
they become only mildly suppressed at the experimentally accessible energies. These opera¬ 
tors would, in particular, modify the dispersion relations of elementary particles producing 
species-dependent shifts of their maximal velocities [28]. Experimental constraints on this ef¬ 
fect are very tight [24] and the logarithmic suppression is by far insufficient to satisfy them. 
Higher-dimension Lorentz violating operators are snppressed by A* and are less problem¬ 
atic. Still, important bounds on their contributions into particles’ dispersion relations can 
be obtained from astrophysics [29, 30, 31, 32, 33]. 

There are two ways to address this challenge. One is to stay within perturbative regime 
and assume that the theory has another symmetry which entails LI as an accidental symmetry 
at low energies. Remarkably, the role of the extra symmetry can be played by non-relativistic 

^Lorentz violating terms of dimension two cannot be constructed from the fields of the theory and those 
of dimension three can be forbidden by imposing the discrete CPT symmetry. 

^Of course, this does not imply that any non-relativistic theory becomes LI in the infrared; for example, 
there are RG flows terminating at fixed points with Lifshitz scaling. 

^In principle, one could consider a setup where at tree level LI is an exact symmetry of the matter sector 
and is violated only in gravity [26]. Then, Lorentz violation is transmitted to matter only through loops 
with gravitons which are suppressed. However, the departure from the homogeneous Lifshitz scaling in the 
UV required by this approach can compromise the high-energy properties of the theory [27]. 
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supersymmetry [34, 35]. Of course, in a realistic model the supersymmetry must be broken 
which leads to the generation of potentially dangerous Lorentz violating operators. However, 
if the scale of supersymmetry breaking is hierarchically lower than A*, the coefficients of these 
operators are strongly suppressed and can be compatible with the experimental bounds. It 
should be mentioned though that it is unclear if the required version of supersymmetry is 
compatible with Lifshitz scaling at high energy [36]. 

Another approach was proposed in [21] (see also [37]) and is illustrated in Fig. 1. It is 
based on the observation that in strongly interacting theories the RG flow (understood in 
the Wilsonian sense) is fast and generally results in a power-law dependence of the couplings 
on the RG scale. Thus, if the Standard Model is embedded into a theory that passes 
through a strongly coupled RG evolution between the Lorentz violation scale A* and some 
lower scale A/^j (still above the experimentally accessible energies), the coefficients of the 
Lorentz violating operators are suppressed by (A//j/A*)“ with a > 0. The latter exponent 
is controlled by the dimension of the Least Irrelevant Lorentz Violating Operator (LILVO) 
[21]. If a and the hierarchy between the scales Ajji, A* are large enough, the experimental 
constraints will be satished. Note that in this scenario it is sufficient that strong coupling 
occurs only in the matter sector and therefore is described, at least in principle, by the 
standard methods of quantum held theory. Gravity can remain weakly coupled at all energies. 

It is worth stressing that this scenario implies that the whole SM sector passes through 
the regime of strong coupling. In other words, all SM helds must emerge as bound states 
of the strong dynamics with the compositeness scale set by A/^j. Gonstruction of realistic 
models with this property is challenging. A useful tool is provided by the gauge/gravity 
(holographic) correspondence that enables to describe 4-dimensional strong dynamics in 
terms of weakly coupled theory in one dimension higher. This approach was used in [21] to 
implement the above mechanism for emergence of LI in a scalar toy model. The purpose of 
this work is to extend the analysis to theories with fermions. We focus on the case when 
the non-relativistic UV theory is invariant under spatial rotations, so fermions are dehned 
as helds transforming in the spinor representation^ of SO{3) (or SO{d — 1) in the general 
case of (d — 1) spatial dimensions). 

Inclusion of fermions in the framework of emergent LI brings the following puzzle. The 
key role in the structure of the Standard Model is played by chiral fermions. However, the 
very dehnition of chirality relies on the existence of the Weyl representations of the Lorentz 
group. How chiral fermions can appear in a theory that fundamentally does not possess 
LI ? We will see that the answer to this question relies on a simple kinematic property: an 

"‘We shall assume that these fields obey anti-commutation relations. 
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Figure 1; Various scales in the proposal for emergence of Lorentz invariance due to strong 
coupling. Above A* the theory is essentially non-relativistic. The strongly coupled RG flow 
below A* drives the theory towards a Lorentz invariant hxed point. A relevant deformation 
terminates the flow at A/r where the theory enters into a conhning phase and the Standard 
Model helds emerge as composite states. 


SO {3) spinor of the UV theory has two independent components which matches precisely 
the number of components of a four-dimensional Weyl spinor^. The degrees of freedom of the 
UV fermion are dressed by the strongly coupled RG evolution and at low energy constitute 
a chiral fermion. A necessary and sufficient requirement for this to happen is presence of 
a gapless mode in the low-energy theory. We will present holographic models of RG flows 
where this is indeed the case. These constructions will enable to compute explicitly the LV 
contributions in the low-energy theory and verify their power-law suppression. 

The paper is organized as follows. In Sec. 2 we briefly review the holographic formalism 
for fermions in the standard relativistic setting. In Sec. 3 we consider a simplihed model 
describing a non-relativistic fermion coupled to a strongly interacting relativistic sector. The 
latter is taken to be vector-like, i.e. it does not contain any massless chiral bound states. We 
show that this construction leads to appearance of a gapless chiral mode in the low-energy 
theory together with emergent LI. In Sec. 4 we turn to a setup modeling an RG flow from 
a theory with Lifshitz scaling in UV towards a LI infrared hxed point. We assume that the 
RG how is stopped by the theory entering into a conhning phase. We show that a suitable 

^The match holds for any even-dimensional space-time. 
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pattern of confinement again leads to emergence of a chiral bound state. In Sec. 5 we discuss 
emergence of LI in theories with fermions in odd space-time dimensions. Section 6 is devoted 
to conclusions. 


2 Holography for relativistic fermions 

In this section we review the holographic correspondence for fermions in the LI case following 
[38, 39, 40, 41]. It establishes a relation between a strongly coupled conformal held theory 
(CFT) with large number of degrees of freedom in d dimensions and gravity in [d + 1)- 
dimensional anti-de Sitter (AdS) space-time. The latter has the metric. 



where I is the AdS radius. A fermionic operator O^p in the CFT corresponds to a fermion ip 
in AdS. The action for a free fermion with mass M reads, 

S' = - y d'^^^x^/\^\ i{iplpip - Mipip) + Sa , (2) 

where ip = ■^^F* and the Dirac operator is® 

U) = e^V^Du , Dm = = ^[F^, F^j ; (3) 

F"^ are the (d -|- l)-dimensional Dirac matrices. The boundary term Sa is needed to ensure 
the correct variational principle and will be specihed shortly. The sign of the action is hxed 
by unitarity [41]. Choosing the diagonal vielbein and computing the corresponding spin 
connection, 

= {u/l) = r]pju , (4) 

with all other components of u vanishing, one simplihes the expression for the Dirac operator, 

?/ n 

lp = j{T^du + T>^d,)--T'^. (5) 

®We use capital Latin letters from the beginning of the alphabet for indices in the local Lorentz 

frame in d -b 1 dimensions; capital letters from the middle of the alphabet M,N,... are used for (d + 1)- 
dimensional space-time indices; Greek letters ... denote indices in d dimensions spanned by t and cc*, 
i = 1,... ,d— 1 (as this space-time is flat, we do not distinguish the tangent-space indices and the indices in 
the local Lorentz frame). 
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To solve the Dirac equation 


( 6 ) 


[Ip — M)'iIj = 0 

one decomposes the spinor into eigenvectors with respect to T^, 

^ = + r>± = ±^±. (7) 

Performing the Fourier transform along the d flat dimensions, ip oc we obtain the 

general solution, 

ip-{p,u) = {pup^+^^/^{xi{p)JMi+i/ 2 {pu) + X 2 {p)Ymi+i/ 2 {pu)) , (8a) 

i^+{p,u) = {pup‘^^^'>^^^-^^^{xi{p)JMi-i/2{pu) + X2{p)Ymi-i/2{pu)) , (8b) 

where we have introduced p = y/—pP^- 

For Ml > 1/2 the mode associated with X 2 is not normalizable. It is interpreted as a 
source for the operator in the CFT. More precisely, the source y is dehned by the relation 

X = \im(u/l)^^-^^^ip. . (9) 

n ->-0 

The partition function, 

Zlx] = j dipdtp , ( 10 ) 

with the integral taken over conhgurations satisfying the boundary condition^ (9) provides 
the generating functional for the correlators of O^. The CFT without sources corresponds 
to X 2 = 0. The dimension of is then read off from the small-u behavior of the component 

ip+ oc dimO^ = Ml + ^ . (11) 

P 2 

At Ml < —1/2 the situation is reversed. Now the non-normalizable mode is associated with 

Xi which dominantly contributes into the component ip^. Thus, the latter is interpreted as 

a source, whereas ip- determines the dimension of the dual operator, 

dimO; =-M/ + ^ . (12) 

In the window —1/2 < Ml < 1/2 both xi and X 2 modes are normalizable and the system 
admits two different quantizations [41]. In this paper we restrict to the range Ml > —1/2 

^Generally, the integration in (10) should be performed over all fields present in the AdS theory, including 
the metric, subject to appropriate boundary conditions. For simplicity, we omit them in our discussion. 
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and choose the quantization with the source related to ip- by Eq. (9); the dimension of 
the operator is given by (11). Note that at Ml = —1/2 the expression (11) saturates the 
lower unitarity bound on the dimension of a spin-half operator in d-dimensional space-time 
[42, 43, 44], 

The variation of the total action (2) must vanish on the solution. Integrating the bulk 
term by parts and setting the contribution at m = cxd to zero we obtain, 

5S = limi [ dfx (e/l) ^ {ip-5ip^ — ip+5ip-) + 5 Sq . (13) 

^^0 Ju=e 

The component tp- is subject to the Dirichlet boundary condition (9), so its variation van¬ 
ishes. However, 5ipj^ cannot be set to zero. To cancel this contribution one chooses [40] 

Sq = — limi [ d^x (e/l) ^ i’-'tp+ ■ (14) 

^^0 Ju=e 

If the AdS space is bounded by a brane at hnite e, which corresponds to a hnite UV cutoff 
in the dual CFT, the term (14) must be evaluated on this UV brane. 

We assume that the CFT enters into a conhning phase at low energy, so that it possesses 
a discrete spectrum of excitations, cf. [45]. This is realized by placing a brane® at u = L 
(referred to as “IR brane” below) and removing the region u > L. Physically, one interprets 
L~^ as the conhning scale. We will consider two choices of boundary conditions on the IR 
brane, which we call and B-. In the case one imposes vanishing of ip+ on the brane, 

B*- (15) 

while the tp- component is arbitrary. In the case B- the situation is opposite. 


B- : tp- = 0 , tp^ — arbitrary . 


( 16 ) 


The spectrum of masses is obtained by substituting (8) into Eq. (15) or (16) (recall that 
X 2 = 0 due to the boundary condition at m —)■ 0), 



.-1 (MZtI/2) 
fi'j , 


(17) 


where is the j-th positive root of the Bessel function. Massless modes should be studied 
separately. A straightforward analysis shows that massless modes are absent in the case B+. 
On the other hand, in the case B- there is a massless mode of the form. 


i’+ = , V'- = 0, 


( 18 ) 


®To solve the Einstein equations, the brane must have negative tension. This does not lead to pathologies 
if the gravitational field obeys suitable boundary conditions, see e.g. [46]. 
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where the spinor xo satishes the equation P/jT^Xo = 0 . 

Henceforth we specify to the case of even d] we will return to the odd-d case in Sec. 5. 
It is convenient to choose the {d + l)-diniensional gamma matrices as 

T>^ = Y, /i = 0,...,d-l , r“ = 7'^+^ (19) 

where 7 ^ are d-dimensional Dirac matrices and 7 ^^+^ is the chirality matrix — the d-dimen- 
sional analog of 7 ^. From (9) we see that the source y and the corresponding operator 
O^p belong to the Weyl representation of the d-dimensional Lorentz group. However, the 
spectrum of the bound states contains a chiral massless mode only in the case i3_. In the 
case all eigenmodes are massive and have both left and right components combining into 
the full Dirac representation. In this sense the dual low-energy theory arising in the case 
is vector-like. 


3 Chiral modes from non-relativistic fermions 

We want to couple the (deformed) CFT of the previous section to a non-relativistic fermion 
and study the resulting RG flow. To this end we introduce a UV cutoff A* and at this scale 
promote the source y to an elementary dynamical held with Lorentz violating action S^. 
The total action at the scale A* takes the form, 

^ = Soft + S^ + A;^' J d^x y Op, . ( 20 ) 

We will assume to be invariant under spatial rotations®. For concreteness, let us start 
from a simple choice, 

Sx = -b j d^xi{xi^doX + vxYdix) , ( 21 ) 

where the velocity v of the y-held is different from unity. The parameter b with dimension 
of length has been introduced to render the action dimensionless. Note that we have not 
included the Lorentz violating operator without derivatives y 7 ®y which is odd under the 
charge-conjugation (and CPT). For the moment we will assume these symmetries, though 
they are not essential for the present setup (see below). 

On the AdS side the above setup is realized by cutting the space-time with a UV brane 
at M = A^^ and supplementing the action ( 2 ) with the boundary term, 

Suv = ~b [ d‘^xi{'4)_x^do'ip_ -1-w'?/'_7*(9j'^_) . (22) 

J 

®It is worth noting, however, that the analysis of this section can be straightforwardly generalized to the 
case without spatial isotropy. 



In this way we obtain a theory in the slice of AdS bounded by a UV and IR branes, cf. [46, 45]. 
Without loss of generality one can identify with the AdS radius /; this is done from 
now on. 

The canonical dimension of the elementary fermion y is (d—1)/2. Recalling the dimension 
(11) of the operator we observe that the interaction between y and in (20) has 
dimension d + Ml — 1/2. If Ml >1/2 this interaction is irrelevant and the Lorentz violating 
sector simply decouples from the CFT in the infrared. Hence, no LI emerges in this case, 
similar to what happens in the scalar model at irrelevant coupling [21]. In what follows we 
restrict to the range 

-1/2 <Ml < 1/2 , (23) 

when the interaction is relevant and generates a strongly coupled RG flow. Furthermore, 
we will focus on the choice of the IR boundary conditions. As discussed in the previous 
section, the CFT in this case does not possess any massless modes. Then the total system 
has a single gapless mode coming from the elementary fermion y which, as we are going to see 
shortly, survives down to IR. On the other hand, in the case of boundary conditions the 
total number of gapless modes is 2 (one from the CFT and one from y). It is straightforward 
to show that they pair and form a gap [45]. As our aim is to study the emergence of chirality, 
this case is not of interest to us. 

From (2), (22) one derives the boundary condition on the UV brane, 

[5(7°(9o + V 1^^^ = 0 . (24) 

The spectrum of eigenmodes is determined by imposing (24) together with the IR boundary 
condition (15) on the general bulk solution (8). This leads to the wavefunctions, 

= Xif-iu) , = ^^^yi/+(M) , (25) 

P 

with 

f±{u) = JmitiMpu) - Ymi^i/2{pu) (26) 

L Ymi-i/2{pL) \ 

and the spinor Xi obeying the relation, 

Xi = 0. (27) 

We are interested in the gapless mode whose d-dimensional momentum squared is much 
smaller than the conhnement scale, 

pi <^pL <^1 . (28) 
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Upon expanding the ratio f-{l)/f\{l) under these assumptions Eq. (27) simplihes, 


a;7°(l + (l-2M/)yQ 


1-2MI 




1-2MI 


Xi = 0, (29) 


where we denoted u = —po, h = Pi- It is natural to take the parameters b and I to be 
of the same order as they both characterize the UV properties of the system. Then (29) 
has the same form as the standard relativistic equation for a Weyl spinor, up to corrections 
suppressed by a power of the small ratio {l/L) of the IR and UV cutoffs. Note that the 
maximal suppression achievable in this setup is (//T)^ which is reached at Ml ^ —1/2. 
Squaring the combination of the y-matrices in (29) yields the dispersion relation, 


u = ±|k 


1 + (n 


1)(1 - 2Ml)j 



(30) 


We see that the velocity of the mode is close to one for large hierarchy between the UV and 
IR scales. The equation for the amplitude of the positive-frequency mode then takes the 
form, 

SkXi = 0 , Sk =-|k|7°/cjyV (31) 


Recalling that xi also satisfies the chirality condition T^y = 7 ‘^’''^yi = —yi (see (7), (19)), we 
conclude that it has independent components with the same structure as a Weyl spinor 

in d dimensions. In particular, in the case d = 4 it has a single component corresponding to 
the spin pointing along the direction of the spatial momentum. 

The above derivation was carried under the assumptions (28). Using (30) one hnds that 
the rightmost inequality is satished as long as the momentum of the mode does not exceed 
certain upper limit. 


k| < 




(32) 


This bound is parametrically larger than 1/L, but smaller than l/l. To understand what 
happens at higher momenta we expand f±{l) for pi 1 treating the product pL as a quantity 
of order one. This yields, 

piUU = )j(pr)(pi)'-“". 

where J7 [z) is a combination of Bessel functions with order-one coefficients. The eigenvalue 
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equation following from (27) can be cast into the form^*^, 


^plY+2Ml 

J{pL) 


2{v 



(33) 


Note that the l.h.s. depends only on the LI combination p = \/a;^ — k^. Inverting (33) we 
obtain the dispersion relation, 


u 


2 








1-2MI 


(34) 


where F{z) is a dimensionless function. The dispersion relation at low momenta is obtained 
by expanding F{z) in Taylor series. For the gapless mode the zeroth-order term vanishes 
and the linear term gives (30). It is instructive to consider the next term in the expansion 
which gives a contribution into the dispersion relation of order, 


(n-1)2(6//) /2k^ 


(35) 


For Ml < 0 the effective scale entering this contribution is parametrically higher than the 
UV scale /“L Thus, we conclude that the strongly coupled RG flow can suppress also the 
higher-order Lorentz violating corrections to the dispersion relation. For high |k| the Taylor 
expansion does not apply. The analysis then depends on the sign of (n — 1). If the velocity 
of the elementary fermion is superluminal, n > 1, Eq. (33) implies that at |k| oo the 
combination pL approaches the hrst positive root of the function J7, 


pLzz zi , J{zi) = 0 . 


In this regime the dispersion relation again becomes relativistic, but with a non-vanishing 
mass, cf. [21]. It is unclear whether this recovery of LI dispersion relation at high momenta is 
a peculiarity of holographic models, or is a more general property of emergent LI in strongly 
coupled systems with superluminal propagation. On the other hand, in the subluminal case 
n < 1 the deviation of the dispersion relation from the relativistic form grows with |k|. Both 
in the sub- and superluminal cases it is easy to check that the spinor wavefunction satisfies 
(31) for all k and hence describes a state with a fixed helicity. 

Let us return to the case of low momenta satisfying (32) and write down the effective 
action for the gapless mode. This is obtained by substituting the wavefunctions (25) into 

deriving (43) we have neglected terms on the l.h.s. suppressed by positive powers of pi. This is 
justified only for the gapless mode which we are interested in. For higher excitations the omitted terms are 
important. 
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the action composed of (2), (14) and (22). One finds that the bulk contribution vanishes 
and the action is given by the boundary terms on the UV brane^^, 


Seff = / 




P 


XiPiiYXi + b\f-il)\ Xi{-(^X +vka^)xi 


= -i d'^x 


hrd^Xi + (1 - 2M/) 


b/I 


l\L 


1-2MI 


Xi(7°c>i) + l'7'9i)Xi 


(36) 


where in the last equality we changed to the canonically normalized held 


Xi = Xi 


P 


(37) 


and switched back to the coordinate representation. We observe that the action consists of 
the LI piece describing a Weyl fermion and a suppressed Lorentz violating contribution^^. 
Variation of this action reproduces the equation of motion (24). 

From the viewpoint of the dual held theory these results can be understood as follows. 
The coupling to the elementary fermion generates the RG how towards the CFT correspond¬ 
ing to the alternative quantization of the bulk theory mentioned in Sec. 2. The latter contains 
a fermionic operator with the dimension (12). The deviation from LI is governed at low 
energies by the irrelevant deformation of the CFT action, 


SScFT OC j d^x . (38) 

In the presence of conhnement this deformation produces Lorentz violating corrections to 
the ehective action of the bound states which depend on the ratio I/L precisely in the same 
way as in (36). 

The analysis can be easily generalized to the case of an arbitrary bare action for the 
elementary fermion. Consider, for example, replacing (21) with 

J + 7°^(-A))x , (39) 

where r2(—A) is a function of the spatial Laplacian A = didi. This action violates the charge 
conjugation and CPT. In the absence of coupling to the CFT it describes a non relativistic 
fermion with dispersion relation 

a; = C(k2) , (40) 

this stage we do not impose the UV boundary condition (24). 

^^For a single field this contribution can be removed by a rescaling of the space or time coordinate. 
However, this will no longer be possible if the system contains several particle species. 
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whose spin degree of freedom is decoupled from the translational motion. In d = 4 this 
is a fermion with two states corresponding to the spin projections ±1/2 on an arbitrary 
axis. However, the situation changes qualitatively once we couple it to the bulk theory. The 
expression (39) translates into a boundary term on the UV brane, 


S'uv = b 


' U=l 


ct’-x iIj-[ - i'y^do + , 


which leads to the effective action for the light mode, 


= -i (fx 


± (1 - 2MI) 


hf I 


l\L 


1-2MI 


{do + 


(41) 


(42) 


The non-relativistic term in the action is again suppressed relative to the LI kinetic term 
induced from the bulk. The corresponding dispersion relation reads. 


u = ±|k| 


1 - (1 - 2MI) 


bf I 


l\L 


1-2MI-I 


+ - 2MI) 


hf I 


l\L 


1-2MI 


(43) 


and the positive-frequency eigenspinor xi satishes the relation (31) implying that its polar¬ 
ization is aligned with the momentum. Note that in this setup CPT emerges^^ together with 
LI. 


4 Emergent chirality from Lifshitz flows 

4.1 The domain wall geometry 

Lifshitz space introduced in [47, 48] is a useful tool for holographic studies of non-relativistic 
theories. It possesses the isometry 

11—)■ Xi^Xxi, i = l,...,d —1 (44) 

u H-)■ Am, (45) 

where A is a scaling parameter. The AdS space-time is a special case of the Lifshitz geometry 
with z = 1. By extrapolating the ideas of holography to other values of one expects the 
physics of Lifshitz space-time to capture the properties of strongly coupled d-dimensional 
non-relativistic theories invariant under the anisotropic scaling of time and space (44). Em¬ 
beddings of Lifshitz solutions into supergravity and string theory have been discussed in 
[49, 50, 51, 52]. 

12(0) is non-zero, there is a small violation of CPT for modes with very low momentum, which can be 
interpreted as the presence of a chemical potential of order 12(0)(Z/L)^“^'^5 
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To describe an RG flow from a Lifshitz UV fixed point to an IR theory with emergent LI, 
we need a geometry interpolating between Lifshitz space at small values of the holographic 
coordinate, m —)■ 0, and AdS at m —)■ cxd. Geometries with these properties can be obtained 
as solutions to the Einstein equations with matter represented by a massive vector held Vm- 
The action reads 

= ~ . («) 


where Ac < 0 is a cosmological constant and Fmn = QmVn — QnVm is the held strength. A 
detailed study of this system was performed in [53] (see also [21]). It was shown that it has a 
family of solutions labeled by a parameter /*, that are static and invariant under SO{d — 1) 
rotations of the spatial coordinates. The metric and the vector held have the form, 

^ dxidxi + g‘^{u)du‘^), Vt = -T^f{u)j{u), Vi = Vu = 0 , (47) 

\uJ Mvu 


where 


,2 d{d — 1) 

2Ac 


( 48 ) 


The functions /, g, j depend on u through the combination u/l^ and can be found numeri¬ 
cally. They have the asymptotics 


/ = fo{u/kY \ 9 = 9o , 3 = 3o , M 0 , (49a) 

/ = 1 + , S' = 1 + 5-00(m//*)"^"'" , j = 3oo{u/h)~'^^ , u^oo . (49b) 

Here /o,oo, fi'o,oo, jo,oo are constants of order one and 

av = -- + {Mylf . (50) 

For a given /* the solution has the form of a domain wall interpolating between the Lifshitz 
and AdS spaces and centered at m = /*. In Fig. 2 we plot the functions j{u) and g{u) 
for the case d = 4 and several values of the Lifshitz exponent. Numerical analysis shows 
that all three functions /, g, j are monotonically decreasing with u, so that the coefficients 
foo, goo, joo in (49b) are positive. 

The domain wall solutions exist for the vector field masses in the range 

( 51 ) 
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g 


j 




Figure 2: Functions j{u) and g{u) describing the domain wall solution interpolating between 
Lifshitz space at u —)■ 0 and AdS at n —)■ cxo for several values of the Lifshitz exponent and 
d = 4. From top to bottom: z = 3, 2.5, 2, 1.5. The position of the domain wall has been 
set to unity. 

The Lifshitz exponent z is related to the space-time dimensionality d and the quantity Myl. 
We will not need the explicit expression, which can be found in [21], and just note that for all 
values of the vector held mass satisfying (51) the exponent lies in the range 1 < z < d — 1. 
Also, due to the upper bound (51) on Myl the exponent ay is bounded from above and 
numerically turns out to be quite small for interesting values of d (for example, for d = 3 
and d = 4 the upper limit on ay is 0.13 and 0.35 respectively). 

In the dual picture the position of the domain wall sets the energy scale A* = where 
the RG how makes transition between the vicinities of the Lifshitz UV theory and the LI 
IR hxed point. The approach to the latter is governed by the vector operator Oy dual to 
the bulk vector held. The dimension of this operator in the IR theory is found using the 
standard rules of AdS/CFT, 

dim Oy = d + ay . (52) 

Thus, the vector operator is irrelevant and the IR hxed point is attractive as long as ay > 0. 
From (50) we see that this condition coincides with the lower bound on the vector held mass 
in (51). 
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4.2 Fermions in the Lifshitz flow 


We now add fermions to the above setup. We consider the action, 

S = - j d’‘+^x^/\i\ ii>{$ - 0"'’ - M)i<, 


(53) 


where the Dirac operator and covariant derivatives are given by (3) and we have introduced 
a direct coupling of the fermion to the vector field via the operator^"^ 

^(V) ^ iM^VNVMe^T^D^ . (54) 

It modifies the effective metric felt by the fermion in non-zero vector held background, 

Qmn '-t Qmn — i ^— VnVm ■ (55) 

The coupling constant ^ can be different for different fermion species and serves to probe 
species-dependent properties of the RG how. Note that we have omitted the boundary term 
Ss- In contrast to Sec. 3, we will concentrate here on the behavior of the dual held theory 
without sources and correspondingly impose vanishing boundary conditions on at m —)■ 0. 
In this case the boundary term vanishes together with its variation. 

Taking a diagonal vielbein and computing the non-vanishing components of the spin 
connections, 

f' f A.. 

(56) 


^liOO — ? ^uij 

9 ug ug 


m 
T) • 


one obtains the expressions for the operators Ip and Ip in the background (47), 


JP = 


Ip^^^ = -f 


u 


du + 


Iff 


2\f 


r“ + + jT% , 


u 




IL _ i 


(57a) 

(57b) 


If 2lg\f u 

Decomposition of the spinor into eigenstates of r“ and Fourier transform along the {t, x 
coordinates yield the system of equations. 


±((9„ -f F±{u))ij± + ig{- a; 7 ° -h ka")'P^ = , 


where 


F+iu) = 


H + Unf d + U^ 


f 


u 


Mlg 


u 


G{u) = g(l- 


1 + er 

/ 


(58) 


(59) 


do not consider a coupling without derivatives ildYNe^T^ip which can be forbidden by the symmetry 
Vn !->■ —Fat. 
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Note that the l.h.s. of (58) contains only LI operators, while all Lorentz violating contribu¬ 
tions have been grouped on the r.h.s. 

To obtain a discrete spectrum of excitations, we cut the space-time by an IR brane^® 
placed at M = L 3> /*• As discussed in Sec. 2, in the relativistic case there is a massless 
fermion mode for the choice of boundary conditions (16). Let us show that this choice 
also gives rise to a gapless mode in the Lifshitz flow geometry. Consider hrst the limit of 
vanishing frequency and momentum. Setting ca = A;* = 0 in (58) we obtain the solution, 

f-L 

= Xoexp / F^{u')du' , = 0 , (60) 

J U 

where Xo is a constant chiral spinor satisfying^® = Xo- This solution behaves as 

oc («/L)“^+'^A 

in the AdS region u ^ /*. Taking into account the AdS measure f duu~^ entering into the 
mode normalization, we see that the normalization integral for (60) is saturated at m ~ L if 
Ml > —1/2. In other words, the mode (60) is suppressed in the Lifshitz part of the space- 
time. This suggests the following strategy to solve Eqs. (58) for non-zero momenta. As the 
zeroth approximation one takes (60) together with the LI dispersion relation^^ a; = |k|. This 
yields the constraint on the spinor amplitude, 

SkXo = 0 , (61) 

where has been dehned in (31). Thus Xo corresponds to a hxed helicity. The Lorentz 
violating contributions are then treated as small corrections. 

We write, 

-|- , ifj- = , ca = |k| -|- . (62) 

Substituting into (58) and keeping only terms linear in the perturbations we obtain, 

{d^ + + tg = 0 , (63a) 

— {du + F_)^jJ^y + i5(Skt/’+^ — = — i|k|G7°'^® . (63b) 

One multiplies the first equation by Sk and uses the identity = 0 to eliminate the second 
term. The remaining equation has the solution, 

Sk^/’+^ = Xexp/ F+{u')du\ (64) 

J U 

^^The brane energy-momentum tensor required for a static solution of Einstein’s equations is composed 
of a negative tension and a contribution satisfying the null energy condition [ 21 ]. 

^®Recall that for our choice of the bulk E-matrices E” = 7 '^“'"^, see (19). 

^^We focus on the positive-frequency mode; for a; = — jkj the analysis is similar. 
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where the constant spinor x satishes the relations 

= -X , SkX = 0 . 

Substitution of (64) into (63b) yields the solution for 

(u) =(u) exp / F_{u')du', 

J U 

where 

y^iu) = [ du' [i(|k|G'(M') - y^y{u'))yxo + i 9 W)x\ exp 


-2MI 


g{u'')du'‘ 


u'‘ 


(65) 


( 66 ) 


(67) 


In deriving the last expression we have used — F_ = —2Mlg/u. Note that the outer 
integral in (67) is taken from m = 0 to ensure vanishing of the held on the UV boundary. 
We now impose the boundary condition B- on the IR brane, 


^_(L) = 0. 

Multiplying (67) successively by Sk and Sk = |k| 7 ° + kiY and using the relations, 
Sk7°Xo = 2|k|xo , Sk7°Xo = 0 , SkX = 2|k|7°x . 
we hnd that the spinor x must vanish, 

X = 0 , 

whereas the correction to the dispersion relation reads. 


( 68 ) 


(69) 


(70) 


y du Giu) exp 

, .( 1 ) _ 117 

1 

to 

1 _ 1 

uJ — 1 Jv| ^ 

Jg du g{u) exp 

Ju u' 

Finally, the correction yy is found from (63a), 


yy = yy iu) exp 

f Fj^{u’)du' 

' U 


(71) 


with 


ijjy{u) = —i I du'g{u')T^k.yy{u')exp 


(1)7.,o 


2Ml 


g{u")du" 


u" 


(72) 


(73) 


Note that the corrections preserve the property that the fermion wavefunction on the IR 
brane 'ip+iL) satishes Skt/’+(7/) = 0 implying that it describes a particle of hxed helicity. 
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Our next task is to infer the scaling of with the IR cutoff L~^. We notice that for 
Ml > —1/2 the integral in the denominator of (71) is saturated in the AdS region where 
g{u) can be replaced by one. This yields, 


du g{u) exp 


-2MI 


' g{u')du' 


u' 


1 + 2MI 


(74) 


To evaluate the numerator we recall that the functions f,g,j depend on u only through the 
combination 2 : = u/U (see Sec. 4.1). Switching to the latter variable inside the integrals we 
have, 


du G{u) exp 


-2MI 


" g{u')du' 


u' 


rzo 


pL/h 

= h dzG{z)exp 

Jo 

r^o 


-2MI 


g{z')dz' 


-2MI 


= I dz G{z) exp 

fL/U 

■ / dz G{z) exp 

J ZQ 

/ j \ 2av — l 


Z' Jzo Z’ 


Z' 

cL/U 


2Ml j 

\ 2Ml 


g{z')dz' 


z' 


where the coefficients 
/oo - iJ 


(X\ — 


2 

CXD 


2Ml - + 1 ’ 

rzQ 

, Ti/f} I 
02 = Zq 


/ dzG{z) exp 


/o 


2Ml 


g{z')dz' 


z' 


(/oo - WDZq 




2Ml — 2ciy T 1 


(75) 

(76a) 

(76b) 


are independent of and L. On the other hand, they depend on the direct coupling ^ 
between the fermion and the vector held and thus can be different for different fermion 
species. In the second line of (75) we divided the integral in two parts with the intermediate 
point Zq lying in the range 1 2:0 T//* and then used the asymptotic form of the metric 

functions (49b) at 2 : > zq. Of course, the expression in the last line of (75) is independent of 
the choice of zq. Combining (75) and (74) we obtain 


u 


( 1 ) 


k 


/, N 2av 

ai(l + 2M/)f^j +a2(l + 2M/) 



1+2MI-I 


(77) 


This gives power-suppressed correction to the fermion velocity as long as ay > 0, Ml > —1/2. 
To sum up, similar to the model of Sec. 3, we have found that the gapless fermion mode 
has an almost LI dispersion relation with the number and structure of independent degrees 
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of freedom matching those of a Weyl fermion in d dimensions (a single helicity component 
for d = 4). Clearly, this mode is described by an approximately LI effective action with the 
fermion field transforming in the chiral representation of the Lorentz group. 

The form of corrections (77) to the relativistic dispersion relation has a natural inter¬ 
pretation in the dual field theory. At the IR fixed point the theory is described by a CFT 
containing a spinor operator and a vector operator Oy with dimensions (11) and (52). 
The RG flow corresponds to deforming the CFT action by an irrelevant Lorentz violating 
perturbation, 

SScFT = j d"a;(ci l^O^y + C2 , (78) 

where ci ,2 are some dimensionless constants. These two terms generate the two contributions 
in (77) after the theory enters into the conhning phase. Note that a double insertion of the 
deformation Oy is required to affect the fermion dispersion relation, as in the unperturbed 
CFT the three-point function O^O^Oy vanishes due to the symmetry of the bulk action 
under Vm ^ -Vm- 

Finally, let us comment on the domain of validity of the calculation leading to (77). A 
necessary requirement is that the correction (72) to the fermion wavefunction is smaller than 
the wavefunction at the zeroth order (60). From the estimates 

'?/’+~ /Xo ~ (79) 

one obtains the condition 

|k| < L-^ min {(L//*)“'^, (L//*)^/ 2 +m«| _ 

Alternatively, the same condition can be derived from the requirement that the correction to 
the gapless mode energy is smaller than the splitting between the energies of the massless 
mode and the next eigenstate of the unperturbed relativistic Dirac operator. The latter 
eigenstate has relativistic dispersion relation with the mass of order 1/L, so the splitting is 
of order (|k|L^)“^. The upper limit (80) on the momentum is parametrically larger than the 
gap 1/L between the gapless mode and the next bound state. It may even seem that the 
perturbative expansion can be pushed to momenta higher than the Lorentz violation scale 
1//* if q;a > 1 and Ml >1/2. However, we expect that the structure of higher order terms 
in the expansion will restrict its domain of validity to |k| 1//*, similar to the case of scalar 

theory [21]. Anyway, this issue is irrelevant for the setup studied in this section: as pointed 
in Sec. 4.1, for the interesting values of space-time dimensionality the exponent ay is smaller 
than 1. 
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5 RG flows with fermions in odd dimensions 


In this section we consider the holographic RG flow in the case when the space-time dimen¬ 
sionality d of the dual held theory is odd. This affects the form of the bulk T-matrices. 
Instead of (19) we choose, 





H = 0,... ,d — 1, 



( 81 ) 


where 7 ^ are the d-dimensional Dirac matrices. Then the analysis of Secs. 3, 4 goes through 
with minor changes. One concludes that the two-brane model with the action ( 22 ) or (41) 
on the UV brane and B+ boundary conditions on the IR brane gives rise to a gapless fermion 
mode with almost LI dispersion relation. The same holds for the case of Lifshitz domain 
wall with boundary conditions in the IR. Of course, there is no notion of chirality in 
odd dimensions and the fermions transform in the Dirac representation of SO{d — 1,1). The 
absence of the gap is protected by spatial parity which forbids a LI mass term for fermions 
in odd dimensions (see e.g. [54]). Indeed, under the reflection of a single axis a fermion in 
odd dimensions transforms as. 


X(a;°, x\ ...) 1 -^ 7^x(a^°, ...) . 


It is straightforward to check that a mass term mxx changes sign under this transformation. 
On the other hand, the original bulk action (2) or (53) respects parity which in {d + 1) even 
dimensions can be defined as 


'lp{x^, x^, x'^,. 


^ T^+^T^^P{x^,-x\x^,...) 




where = iTOT^... 

A new feature that appears in the case of odd d is a possible mismatch between the 
number of degrees of freedom of non-relativistic and relativistic fermions. Namely, the 
spinor representation of the group SO{d — 1) of spatial rotations can be decomposed into 
the left and right parts, each having dimension 2 *^*^“^^/^, which is thus the minimal number 
of components of a non-relativistic spinor. This does not fit into the Dirac representation of 
SO{d — 1,1) which has 2 ^'^“^^/^ components. Therefore, one does not expect emergence of 
LI in systems with minimal non-relativistic fermions. Let us analyze this question explicitly. 
We take d = 3 and consider the holographic setup with two branes similar to that studied in 
Sec. 3. First, we describe the model from the 3-dimensional viewpoint. The spinor operator 

of the dual CFT can be decomposed into the upper and lower components. 


Oy = 



( 82 ) 
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Taking the 7-matrices in the form, 


7 ° = ias , 7 ^ = (Ji , 7 ^ = 0-2 , (83) 

where cij are the Pauli matrices, we observe that Ou and Od do not mix under purely spatial 
rotations. Hence, it is consistent with spatial isotropy to couple only one of them, say Ou, 
to an external dynamical held. This leads to the action, 

(fxxlfOu + b j d^xXu{ido-VL{-ddj)xu, (84) 

where xu is a single-component minimal SO{2) spinor and r2(—A) is an arbitrary function 
of the spatial Laplacian. Note that this action breaks the discrete symmetries C, P, T and 
CPT as all of them interchange Ou and Ou- 

The holographic description of this model is obtained by making only the upper compo¬ 
nent of the bulk spinor ip- on the UV brane dynamical with the action, 

Suv = bf d^x ip^_ ^{ido - il{-A))ip-^u , (85) 

Ju=l 

whereas imposing the Dirichlet boundary condition on the lower component, 

V'-,dL., = 0. (86) 

From (85), (2) we obtain the second boundary condition on the UV brane, 

[b{ido - n{-A))ip-^u + ^+,c/] = 0 . (87) 

Combining with boundary conditions on the IR brane and substituting the bulk solution 
(25), (26) into ( 86 ), (87) yields, 

[mb{u - +uf+{l)]xi,u = 0 . (88) 

Expansion of the functions f±{l) at small argument gives the dispersion relation of the gapless 
mode, 

u = Q{k^){l-2Ml)-(-\ . (89) 

Clearly, this is not relativistic. In the infrared limit l/L ^ 0 the dispersion relation becomes 
degenerate, a; = 0, implying that, the excitation of this mode does not cost any energy. In 
other words, the system contains inhnite number of states with zero energy different from 
the vacuum. This situation appears rather pathological. It is excluded if we impose from 
the beginning one of the symmetries C, P, T or CPT. 
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6 Conclusions 


In this paper we have studied emergence of LI in strongly coupled held theories with fermions. 
We considered two models. In the hrst model a strongly interacting relativistic CFT is 
coupled at a high energy scale to an elementary fermion with non-relativistic action. This 
generates an RG how exhibiting LI in the IR. The second setup describes a held theory 
howing from a UV hxed point with Lifshitz scaling to a LI IR hxed point. In both cases we 
assumed that the theories enter into a conhnement phase at low energies, which gives rise to 
a discrete spectrum of excitations. We used the holographic duality to describe the strong 
coupling. 

In the main part of the paper we focused on held theories living in even-dimensional space- 
time and analyzed how the notion of chirality arises together with emergent LI. We have 
found that if the theory possesses low-lying fermionic excitations, the latter are described 
by approximately LI ehective actions with the fermion wavefunction transforming in the 
Weyl representation of the Lorentz group. The Lorentz violating corrections to the ehective 
action are power-law suppressed by the ratio between the IR conhnement scale and the 
UV scale of Lorentz violation. The exponent in the power-law is related to the dimension 
of the Least Irrelevant Lorentz Violating Operator (LILVO) in the IR theory, in complete 
analogy with the scalar case studied in [21]. A kinematic property that appears essential 
for such behavior is the coincidence between the number of independent components of the 
non-relativistic fermion operator in the UV theory and the relativistic Weyl spinor. 

In the last section we have addressed the case of odd space-time dimensionality. We 
observed that the setups preserving the kinematic match mentioned above still lead to low- 
energy theories with approximately relativistic fermions which are protected from acquiring 
a mass by spatial parity. We also considered a model that violates the above match and 
showed that in this case the low-lying fermion mode does not exhibit LI. We pointed out 
that this situation is excluded if the theory respects one of the discrete symmetries G, P, T 
or CPT. This is one more manifestation of the deep role played by these symmetries in the 
phenomenon of emergent LI. 

From the phenomenological perspective, our work represents a step towards implementing 
the idea of emergent LI in a realistic particle physics setting. Of course, an important open 
issue on this way is inclusion of gauge fields. As discussed in [21], the holographic description 
of strongly coupled RG flows with gauge fields will presumably require considering a more 
complicated bulk sector containing an interaction of the gauge helds with a dilaton. More 
generally, one can ask whether the gauge symmetry can be an emergent property appearing 
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together with emergent LI. 

Another interesting question concerns the minimal value of LV couplings predicted in this 
framework. We have demonstrated that these couplings are suppressed by a power-law factor 
(A/^/A*)", where the ratio characterizes the “duration” of the RG flow leading to 

the emergence of LI. Clearly, this duration cannot be arbitrarily large. In the models inspired 
by quantum gravity the natural value for A* is at or below the Planck scale. In particular, 
in the context of the Horava’s proposal. A* should be less than 10^^ GeV [55, 6]. On the 
other hand, the scale Kjr is bounded below by particle physics experiments. As mentioned 
in the Introduction, A//j sets the scale of compositeness for the SM helds. Direct searches 
for excited fermionic states^® put a lower bound on Kir at the level of a few TeV [56, 57]. 
Taking for the estimate Km > 10 TeV, we find that the duration of the RG flow is bounded 
from below by Km/K^ > 10“^^. In all holographic examples considered in Ref. [21] and this 
paper the exponent a in the suppression has been found to be smaller than 2. Thus, the 
minimal size of LV operators predicted by these models is of order 10“^^, which is marginally 
compatible with the existing bounds [24]. In this light it will be important to understand 
whether the condition a < 2 is universal or there exist RG flows with emergent LI that 
avoid it. 

To conclude, the scenario of emergent LI provides an interesting interplay between high- 
energy particle physics and precision tests of relativity. An improvement of the experimental 
lower limits on the compositeness scale of SM helds and / or tightening of the constraints 
on LV parameters will be able to test this idea or rule it out. 
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